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The paraconductivity, originating from critical superconducting order-parameter fluctuations in 
the vicinity of the critical temperature in a layered superconductor is calculated in the frame of 
the self-consistent Hartree approximation, for an arbitrarily strong electric field and zero magnetic 
field. The paraconductivity diverges less steep towards the critical temperature in the Hartree 
approximation than in the Gaussian one and it shows a distinctly enhanced variation with the 
electric field. Our results indicate that high electric fields can be effectively used to suppress order- 
parameter fluctuations in high-temperature superconductors. 



I. INTRODUCTION 

Due to their high critical temperature, small coher- 
ence length, and quasi-two-dimensional nature, the high- 
temperature superconductors (HTSC) show a much more 
pronounced and therefore experimentally accessible ef- 
fect of thermodynamic fluctuations in the critical region 
of the normal-superconducting transition. In general, an 
enhancement of the conductivity, denoted paraconduc- 
tivity, is observed in HTSC above T c due to the presence 
of superconducting fluctuations. Outside the critical re- 
gion, in the absence of the magnetic field and for small 
electric fields, the paraconductivity can be explained 
in terms of the Aslamazov-Larkini theory of noninter- 
acting, Gaussian fluctuations. The initial expressions 
for the paraconductivity have been extended for two- 
dimensional layered superconductors, a situation very 
much resembling the crystal structure in the cuprates, by 
Lawrence and Doniach. 2 However, it was shown that the 
fluctuation conductivity may be calculated in the linear- 
response approximation only for sufficiently weak fields, 
when they do not perturb the fluctuation spectrum^ 
Reasonably high values of the electric field E can accel- 
erate the fluctuating paired carriers to the depairing cur- 
rent, and thus, suppress the lifetime of the fluctuations, 
leading to deviation from Ohm's law. In connection with 
the low-temperature superconductors, the nonlinearity 
has been studied theoretically for the isotropic case^ 
and also proven experimentally on thin aluminum films. 6 
The issue of the non-ohmic fluctuation conductivity for a 
clean layered superconductor in an arbitrary electric field 
has been addressed by Varlamov and Reggiani^ start- 
ing from a microscopic approach of Gor'kov 8 for dirty 
isotropic superconductors. Essentially the same depen- 
dence on temperature and electric field has been recently- 
derived, together with generalizations for the case of ar- 
bitrary dimension, based on the analytical derivation of 
the velocity distribution resulting from the Boltzmann 
equation for the fluctuating Cooper pairs. 

The above-mentioned theories do not consider the in- 
teractions between fluctuations, so that fluctuations can 
be described by the Gaussian approximation. Thus, the 
quartic term in the Ginzburg-Landau free energy is ne- 
glected. This approximation is known to hold for tem- 
perature values not too close to the mean-field transi- 



tion temperature, but it breaks down in the critical re- 
gion, since the nonlinear character of the TDGL equa- 
tion cannot be neglected for high densities of fluctua- 
tion Cooper pairs. Several worko 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 
have included the interaction between superconducting 
fluctuations in the critical transition region within dif- 
ferent theoretical approaches. The simplest and most 
used one is the Hartree approximation which treats self- 
consistcntly the quartic term in the Ginzburg-Landau 
free-energy expansion. In this way expressions for the 
specific heat have been derived for bulk 10 and layered 11 
superconductors under magnetic field, based on the func- 
tional integral approach. In the frame of the time- 
dependent Ginzburg-Landau (TDGL) theory, Ullah and 
Dorsej*i^ computed the Nernst effect, the thermopower, 
the longitudinal and the Hall conductivity in the linear- 
response approximation for a layered superconductor in 
a magnetic field. Using the same relaxational dynam- 
ics of the TDGL approach, Dorsejsi^ provided expres- 
sions for the fluctuation conductivity in both ohmic and 
non-ohmic regime, for isotropic superconductors of ar- 
bitrary dimensionality and in the absence of a mag- 
netic field. More recentljsi^ the effects of critical su- 
perconducting fluctuations on the scaling of the linear 
ac conductivity, <j(u>), of a bulk superconductor slightly 
above T c in zero applied magnetic field have been in- 
vestigated based on the dynamic renormalization-group 
method applied to the relaxational TDGL model of 
superconductivity, verifying explicitly the scaling hy- 
pothesis cr(w,£) proposed originally by Fisher, Fisher, 
and Huse^ The essential features of the scaling and 
renormalization-group method has been also reviewed 
recently by Larkin and Varlamov^ A more sophisti- 
cated approach based on renormalization procedure and 
diagrammatic techniques has been developed for treat- 
ing the non-Gaussian superconducting fluctuations by 
Ikeda, Ohmi and Tsuneto^ and applied to the longi- 
tudinal conductivity, the magnetization^ and the Hall 
conductivity*^ in the linear-response approximation, for 
a layered superconductor under magnetic field. 

In this paper we shall address the problem of the non- 
ohmic behavior of the non-Gaussian fluctuation conduc- 
tivity for a layered superconductor, a topic that, to our 
present knowledge, has not yet been treated in the lit- 
erature. While the non-linear conductivity under ar- 
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bitrarily strong electric field was derived for a layered 
system^ for Gaussian, noninteracting fluctuations, the 
effect of the critical, strongly interacting fluctuation on 
the ohmic and non-ohmic conductivity was investigated!* 
only for isotropic systems of arbitrary dimensionality, but 
not for the layered Lawrence-Doniach model. The lat- 
ter, however, would be required for comparison to ex- 
perimental data on, e.g., YBa2Cu306+x- The paper is 
organized as follows. In Section [H] the TDGL equations 
are deduced for a layered superconducting system, with 
the explicit consideration of an arbitrarily strong elec- 
tric field, oriented parallel to the layers. The fluctuation 
interaction term is also included in this model, within 
the self-consistent Hartree approximation. Section Ifffl 
presents the resulting equation solutions, obtained with 
the aid of the Green function technique (detailed in the 
Appendix). By considering in detail the necessary cor- 
rection through the UV cut-off procedure, expressions for 
the fluctuation conductivity and the self-consistent equa- 
tion for the renormalized reduced temperature parameter 
are provided. Further, Section HVI treats the limit cases 
of the linear response approximation, the no-cut-off limit, 
and also the isotropic 2D and 3D cases, recovering thus 
results of previous theories. In Section [V] the applica- 
tion of the model is illustrated by comparing the para- 
conductivity obtained in the present theory with that in 
the Gaussian fluctuation approximation for various ap- 
plied electric fields. A comparison between the fluctua- 
tion suppression effects of the electric and magnetic fields 
is also illustrated. Finally, in Section I VII we summarize 
the main conclusions emerging from our analysis. 



II. TDGL EQUATION FOR ARBITRARY 
ELECTRIC FIELD 



For our purpose, we shall adopt the TDGL framework, 
and treat the quartic term in the free-energy expansion 
within the simple self-consistent Hartree approximation. 
The starting point will be the Lawrence-Doniach expres- 
sion of the Ginzburg-Landau (GL) free energy for a sys- 
tem of superconducting planes separated by a distance 
s, with a Josephson coupling between the planes, in the 
absence of magnetic field, 
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where m and m c are effective Cooper pair masses in 
the a6-plane and along the c-axis, respectively. The GL 
potential a — a^e is parameterized by ao = K 2 /2ra^ 
and e = In (T/T ) w (T - T ) /T , with T being the 
mean-field transition temperature and £o the in-plane 
GL-coherence length, extrapolated at T = 0. The critical 
dynamics of the complex superconducting order param- 
eter ip n in the n-th plane will be described by the gauge- 



invariant relaxational time-dependent Ginzburg-Landau 
equation 
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where the pair electric charge is eo = 2e, and the order 
parameter relaxation rate T , given bji22i2ii2£ 
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is related to the life-time of metastable Cooper pairs*!*^ 
T (BCS) _ n h/l6k B (T - T ) through the relation T^ 1 = 
2ar^ BCS \ As it can be noticed, we define the relax- 
ation rate T^ 1 as depending on the actual tempera- 
ture T, while in many other works^t^S^ it is defined 
as a function of the mean-field critical temperature, T Q . 
Of course the difference is negligible near the transition 
point, but one must recall that in the derivation of the 
GL equations from the BCS theory^SLSiS the tempera- 
ture is usually approximated with the critical one, while 
the reduced temperature ]h(T/Tq) is approximated with 
(T - T ) /T . The same happens while the TDGL equa- 
tion is obtained for temperatures near To, as for instance 
in fiefs. |2^and|2^. In some more recent derivations of the 
TDGLj22i2£ however, the original appearance of these pa- 
rameters is preserved, so that the order parameter relax- 
ation rate is written as depending on the actual temper- 
ature. In this mathematically somewhat stricter sense, 
one could therefore write the life-time of the metastable 
Cooper pairs as T iBCS ^> = nh/l6k B T\n(T/T ). 

The Langevin forces £ n (x, t) introduced in Eq. (J2J in 
order to model the thermodynamical fluctuations must 
satisfy the fluctuation-dissipation theorem, and ensure 
that the system relaxes to the proper equilibrium dis- 
tribution. This requirement is fulfilled if the Langevin 
forces £ n (x, t) are correlated by the Gaussian white-noise 
law 

(Cn(x,t)C (x',t')> =2r - 1 fc s T<5(x-x')<5(t-O— , 

' (4) 

where <5(x — x') is the 2-dimensional delta- function con- 
cerning the in-plane coordinates. Since we are interested 
in finding the conductivity for an arbitrary electric field, 
we cannot use the linear-response approximation, so we 
have to explicitly include the electric field in the model. 
In order to compute the in-plane fluctuation conductiv- 
ity, we shall assume the field E along the x-axis (where 
x and y are the in-plane coordinates), generated by the 
scalar potential ip — —Ex. In the chosen gauge, the cur- 
rent density operator along the x direction in the n-th 
plane will be given by 
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so that after averaging with respect to the noise, 
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The time-dependent GL equation J5J) writes 
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dt 
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As already mentioned, the quartic term in the thermo- 
dynamical potential will be treated in the Hartree ap- 
proximation, in the same sense as applied also in pre- 
vious works fi^*22iS£ namely by replacing the cubic term 

^IV'nl VVi m Eq. J7J) with b ( \^p n \ 2 ) ip n . In this way, 



the non-linearity is decoupled, resulting in a linear prob- 
lem with a modified (renormalized) GL potential a — 

a + b (\ip n \ 2 ^, which implies a renormalized reduced tem- 
perature 



e = e ■ 
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The average (\ip n \ J is to be determined, in principle, 

self-consistently together with the parameter e. 

In order to simplify the following computations, we 
shall introduce the Fourier transform with respect to 
the two in-plane coordinates and the layer index, respec- 
tively, through the relations: 
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so that Eq. (JJJ becomes 
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We have int roduced the anisotropy parameter 7 = 
£oc/£o = \J m/m c , with £oc the out-of-plane GL- 
coherence length extrapolated at T = 0. The term 
(k, i) in Eq. I|10[l is the Fourier transform according to 
the rules © of the noise function ( n (x, t) . One can di- 
rectly verify that the following correlation function holds: 



(C, (k, t) Q (k', t')) = 2r - 1 fc B T(2 7 r) 3 (11) 
■5{k~k')6(q~q')S{t-t'). 



Equation (|10|) may be solved using the Green function 
method, as done also by Tucker and Halperin^ and 
Dorsej*i^ in order to solve the relaxational TDGL equa- 
tion with a Langevin noise for isotropic superconductors. 
Our derivation differs however by the fact that it is ap- 
plied for a layered superconductor, and also through the 
differently chosen potential gauge. We denote thus with 
i? g (k, t; k' x ,t') the Green function for the Eq. fTT)|. which 
satisfies 



_ x 8 epT^E d 
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where we have introduced the notation 
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Thus, the solution of Eq. 1)10(1 will then be written 



ip q (\s.,t) = / dt' / dk' x R q (k,t;k' x ,t')C q (k' x ,k y ,t') 



(14) 

The Green function R q (k,t; k' x , t') is computed in the 
Appendix, and given by Eq. (|A.7J) . so that the rela- 
tion l|14|) for the Fourier transformed order parameter 
will write accordingly 



ip q (k,t) = r exp 
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One can notice that the solution 1(15(1 fulfills causality, due to the retarded character of the Green function 
R q (k,t; k' x ,t'). It can be also simplified further, 
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Now, following (jSJ, the current density can be written as 
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while the averaged density of fluctuating Coooper pairs writes analogously 
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where we have taken into account the expression l|ll(l for 
the Fourier transformed noise correlation function, and 
also replaced the notation ai with its value l|13fl . 

Before proceeding further and solve the integrals over 
momentum variables in Eqs. (|17|l and 1)18(1 . we must 
recall the inherent ultraviolet (UV) divergence of the 
Ginzburg-Landau theory, which is not valid on length 
scales less than the zero-temperature coherence length 
£o- The short wavelength fluctuations break down the 
"slow variation condition" for the superconducting order 
parameter, a central hypothesis of the GL approach, 30 
This difficulty can be solved by applying an UV cut-off 
to the fluctuation spectrum, a procedure introduced from 
the beginning in the GL theory^ 



III. SOLUTIONS WITH CONSIDERATION OF 
THE UV CUT-OFF 

The classical^SSi 3 ^ procedure is to suppress the short 
wavelength fluctuating modes through the momentum 
cut-off condition 



k 2 < c£ - 



(19) 



where the dimensionless cut-off factor c is close to unity. 
Also a total energy cut-off was suggested ) 33 i 34 i 35 which 
eliminates the most energetic fluctuations and not only 
those with short wavelengths, 



k 2 + r 2 (£)<<4- 



(20) 



Heuristically, the adequacy of the total-energy cut-off 
can be justified on the grounds of the Gaussian GL ap- 
proach by taking into account that the probability of 



each fluctuating mode is controlled by its total energy 
fi 2 k 2 /2m + age, and not only by its momentumiSiiii Very 
recentlj^ it was suggested that the physical meaning of 
the "total energy" cut-off follows from the uncertainty 
principle, which impose a limit to the confinement of 
the superconducting wave function. It must be however 
mentioned that in the critical fluctuation region, the two 
cut-off prescriptions almost coincide quantitatively, due 
to the low reduced-temperature value e with respect to 
the factor c. 

For simplicity, in the following we shall apply the cut- 
off procedure in its classical form (|19fl on the k-plane 
momentum integrals in Eqs. (|17|l and l]18p. If necessary, 
one can get easily also the result that would correspond 
to the "total energy" cut-off (|20|) . by simply replacing 
c — > c— e. Although for the Lawrence-Doniach model, the 
z-axis momentum q also contributes to the fluctuation 
mode energy with the term {h 2 ^ 2 /ms 2 ) (1 — cosgs), the 
inclusion of a momentum cut-off in this direction is not 
necessary, since the z-axis spectrum is already modulated 
through — tt/s < q < tt/s. One can thus perform the 
integral over g-momentum in Eqs. I|17fl an d (|18|) and get 
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where we used the identity 7 (a;) = (l/2vr) /* r e" os9 J 
for the modified Bessel function Iq(x). 

We shall therefore apply the cut-off by performing the 
momentum integral from Eq. 1(17(1 for k 2 < c£q~ 2 , 
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where we introduced the new dimensionless variable iu = 
H 2 k 2 /2mao and used the first order modified Bessel func- 
tion J a (x) = (l/2n) J^dip cos 93 e 21 ™ 3 ^. The current 
density (|17|l will write eventually, after considering Eqs. 
(|21|l and l|22[) . and introducing the new integration vari- 
able u — 2aoToT, 
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where we have also expressed the parameters To and cto 
with the aid of the in-plane coherence length £o- 



In order to infer the fluctuation conductivity a = j/E, 
we shall use the fact that the function I\ (x) consists only 
of odd argument powers, and satisfies namely the identity 
Ii(x) = (x/2) [Iq(x) — h{x)], so that we can finally write 
the paraconductivity under the form 
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where we have introduced the notations 
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for the anisotropy parameter r, and the characteristic 
electric field Eq, respectively. Relations lll'ol) can be 
expressed also as depending on microscopical parame- 
ters, like the Fermi velocity vp and the electronic inter- 
layer hopping energy J. By identifying the in-plane GL- 
coherence length £o from the microscopic derivation of 
the GL equation in the two-dimensional case, one has in 
the clean limftia £ = (7C(3)/32) 1/2 hv F /Trk B T , which 
implies for T w Tq 



E = 64. 
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where we have also used in the same approximations 
(clean limit and T ~ Tq), the expression for the 
anisotropy parameter r from Ref. l37l as a function of 
the inter layer electron hopping energy J. 

In an analogous manner as presented above for the 
current-density, one can apply the cut-off procedure on 
the momentum integral in Eq. (|18|l . and obtain 
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(26) so that the averaged density of Cooper pairs l|18[> becomes 
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where we passed from r to the variable u, and used the notations One can easily prove that without the UV 

cut-off (i.e. for c — * oo), the expression would be divergent. Replacing c — > c — e would in turn correspond to the 
energy cut-off under the form 1)20(1 . The self-consistent equation (JHJ for the parameter s will write therefore: 



where the factor 



one obtains for the linear response conductivity 



Trh 2 s 



(30) 



has been computed by taking into account the expression 
of the quartic term coefficient b = /i K 2 egfi- 2 /2m 2 , with 
k being the Ginzburg-Landau parameter k = Ao/£o- 

The above Eq. ((24(1 for the paraconductivity, together 
with the self-consistent Eq. (|29|) . both valid for an ar- 
bitrary strong electric field and with the explicit inclu- 
sion of the UV cut-off, are the main results of this pa- 
per. A short comment would be useful regarding the 
application of the cut-off procedure to the k-integrals in 
Eqs. ((27(1 and 1(22(1 . It could seem more appealing to 
apply the cut-off on the translated wave-vector magni- 
tude, such as ky + [k x + (eoET/2h)] 2 < c£q~ 2 , as it was 
approximately done, although in a different gauge, by 
Kajimura and Mikoshiba, 38 while studying the paracon- 
ductivity in arbitrary electric field in the two-dimensional 
case. This would certainly simplify the calculations and 
consequently the factors that depend on the cut-off pa- 
rameter c, but would not really correspond to the actual 
meaning of the UV cut-off, which is to assure the slow 
variation condition for the order parameter by eliminat- 
ing the rapidly oscillating modes in the Fourier expansion 
@. It can be seen for instance that when the dummy 
variable t grows towards oo, the component k x would 
have to approach — oo in order to preserve the cut-off 
condition in the translated form. It can be moreover 
verified that the "translated" cut-off would not give the 
correct result in the E — limit, namely the third term 
in the Eq. ((31(1 below would be missing. 



IV. LIMIT CASES 



A. Linear response limit 
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y/e (e + r) y/[e + c) (e + c + r) 
c(c + e + r/2) 



[(c + e + r)(c + e)] 



3/2 



(31) 



which is the Lawrence-Doniach 2 " formula for the fluctua- 
tion conductivity of a layered superconductor, with the 
inclusion of the UV cut-off. A formally identical expres- 
sion was also inferred by Carballeira et a/., 35 for Gaus- 
sian fluctuations (i.e. with e = e), in order to fit the 
HTSC paraconductivity in the high reduced-temperature 
region. 

The integral in Eq. (gSJl can also be easily performed 
in the limit E — 0, and it yields the relation 
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that matches the formula found by Mishonov and 
Penev?2& with the only difference that in Ref. the 
temperature T is approximated with To in the logarithm 
prefactor. It can be easily verified that Eq. 1(32(1 is also 
coincident with the analogous equation found previously 
by Ullah and Dorsejii 3 . (UD) for the renormalized reduced 
temperature eh, if one considers the zero magnetic field 
limit in the UD formula. In this limit, the sum over Lan- 
dau levels transforms, with the aid of Euler-MacLaurin 
summation formula, to the integral from Eq. 1(32(1 , if one 
takes also into account that the cut-off used by UD cor- 
responds with c = 2. 

We have thus shown that from our Eqs. (O and (HJl 
one can infer in the limit E — -> the already known results 
for conductivity and renormalized reduced temperature 
in the ohmic approximation. 



It is worth comparing to previous results what become 
Eqs. iJUJl and igSJ in the zero-field limit, E = 0. Taking 
into account that /o(0) = 1 and ^(0) = 0, and using the 
integral form for the modified Bessel function Iq (to/2), 



B. No-cut-off limit 

The cut-off procedure is crucial for calculating the av- 
eraged fluctuating Cooper pairs density, since Eq. I(18|) 
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yields a divergent r-integral when one performs the k- 
momentum integrals on the entire k-planc. 

The result (|24JI for the paraconductivity remains how- 
ever finite even if one removes the cut-off (i.e. for 
c — > oo), although it will give then a larger paracon- 
ductivity than in the cut-off case, especially for higher 
reduced temperatures (e > 0.1). Turning back to Eq. 
(|17|l with the cut-off removed, one is allowed to translate 
the integral variable k x so that k x + (cqEt^H) — > k x , 
and after performing the Poisson k-integrals, one obtains 
eventually for the fluctuation in-plane conductivity in the 
presence of an arbitrary electric field E, without consid- 
ering the UV cut-off, the simpler relation 




(33) 

The expression Q33JI is not new. In its form, it is 
essentially similar to the ones found by Varlamov and 
Reggiani and Mishonov et alX for the case of Gaussian 
fluctuations, if one neglects the difference which consists 
in the presence of the renormalized parameter e in Eq. 
(|33|l instead of the reduced temperature e = ln(T/To). 
It must be stated however that Ref. |jj based on a mi- 
croscopical approach^ defines an out-of-plane coherence 
length larger by a factor %/2 than the commonly used 



one. Nevertheless, its result corresponds to our Eq. i|33|) 
when expressed through Eqs. i|2ti|) in the microscopical 
parameters vf and J (the latter is denoted in Ref. by 
w). Mishonov et al~ signalize differences in one of their 
intermediary results with respect to Refs. and 0, but 
they find the same dependence of the non-ohmic conduc- 
tivity on the reduced temperature s and the electric field 
E, by solving the Boltzmann equation for the velocity 
distribution of the fluctuating Cooper pairs. Their for- 
mula differs mathematically from Eq. i.Y.'A) only by the 
presence of two extra-factors before the integral, namely 
the ratio between a- and b- coherence lengths (since Ref. 
Intakes into account also a general in-plane anisotropy), 
and the ratio T/Tq, which in fact comes artificially only 
if one writes the relaxation rate Tq with To instead of T 
(see the comment after Eq. [3} ■ 



C. Isotropic limit 

Results analogous to Eqs. (|24[1 and (|29[1 for the 

isotropic two dimensional (2D) case can be easily derived 
by simply taking the limit £oc = (or r — > 0) and identi- 
fying the interlayer distance s with the film thickness d. 
We will have thus 



(e,E) 



du- 2 



e — In ■ 



T 



IQhd j Q 
irh 2 d 



u e v °' 



dw '< 



/„ ( 2V3~|V 

£/ 



-h 



2^3— u 2 

Eq 



du e 



2\/3— uVr 
Eq 



, (34) 
(35) 



Equations 134|) and l|35|l differ from the analogous re- 
sults of Ref. 1381 because there only an approximate form 
of the cut-off procedure was applied, as we have already 
mentioned at the end of Section IIIII 

If one neglects the cut-off (c — » 00), Eq. (|35|l becomes 
divergent, while Eq. I|34|) takes the already known^2ii4 
form 



T (2D) 
7 NoCut 



(e,E) 



Whd 



du e 



(36) 



with the specification that in Eq. I|36(l the renormalized 
parameter e is present, instead of the reduced tempera- 
ture e. 

If, on the contrary, one preserves the cut-off but takes 
the linear response limit E — > 0, Eqs. I|34|) and l|35|) 
become 



T (2D) 



(?) 



In- 



E=0 

T 



1QM 



1 



1 



2^ K 2 e 2 ^fc B T 1 e + c 
In ■ 



Ti 







irh 2 d 



,(37) 
(38) 



A formally identical expression to Eq. (|37|l is also to 
be found in Ref. |3j| for Gaussian fluctuations (i.e. with 
e — e), while Eq. H38fl i s implicitly contained in the 
results of Ref. HI 



Results for the isotropic three-dimensional (3D) case 
cannot be obtained by just imposing the 3D condition 
s — > (or r — > 00) to Eqs. and because 

these equations were calculated by assuming that a cut- 
off in the z direction is not necessary for layered super- 
conductors. However, in the 3D case, a cut-off for the 
fc 2 -momentum is as necessary as for the k x and k y com- 
ponents. The calculations can be performed according to 
the same scheme as for the layered case, and the results 
are presented here for completeness: 
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r (3D) 



(e,E) 



In 



Tn 



2 poo / 

e Z I 2 -eu-4( 

auu e v 



87T^ Jo 



(39) 



cosh ( 2V3^-u 2 y/w ) sinh h^/Z^-u 2 ^ 



2y/3^-u*\/w 



T 2 f i K 2 e^ k B T 



TT 2 H 2 



du e 



-e u— 4 



sinh 



(40) 



Similarly to the 2D case, if one neglects the cut-off 
(c — ► oo), the r.h.s term in Eq. (|40|) becomes divergent, 
while Eq. 139(1 takes the expression 



T (3D) 
r NoCut 



(e,E) 



_ , dit — — e 

32 % /7r^o Jo Vu 



( 41 ) 

already known£*2*ii for Gaussian fluctuations (i.e. with 
e — e). 

In the linear response limit (i? — ► 0) but with the cut- 
off preserved, Eqs. (|39|l and (|40l) become 



r (3D) 



B=0 



In 



48tt^ 
3 arctan ^y/cfej 



T ^ n 2 e 2 ^k B T 



Ze^/c 5c 3 / 2 



(e + c) 2 (e + c) 2 



7T 2 h 2 



' e arctan 



Equation (|42jl matches thus formally the expression 
found 3 ^ for Gaussian fluctuations (e = e) . 



V. RESULTS OF THE MODEL 

The rcnormalization procedure required for our present 
results consists thus in using the reduced temperature 
parameter e, renormalized by solving Eq. I(29|l . in the 
conductivity expression l|24|) . This procedure causes the 
critical temperature to shift towards lower temperatures. 
In analogy with the Gaussian fluctuation case, we shall 
adopt as definition for the critical temperature T C {E) the 
vanishing of the reduced temperature, e = 0. Thus, in 
the absence of the electric field we can use Eq. Il,i2l) taken 
at T = T c (0) and e = 0, so that one gets 



one can use Eqs. I|29() for any temperature T and field 
E in order to find the actual renormalized e(T, E) , and 
further the conductivity a(T, E). 

In order to illustrate the main features of the Hartree 
approximation for the critical fluctuation model, we shall 
take as example a common material, like the optimally 
doped YBa2Cu3 06+ K - Typical values for the character- 
istic parameters are then: s — 1.17 nm for the interlayer 
distance, £o = 1.2 nm and £oc = 0.14 nm for the zero- 
temperature in-plane and, respectively, out-of-plane co- 
herence lengths, k = 70 for the Ginzburg-Landau param- 
eter and T c q — 92 K for the critical temperature under 
(42V) er y small electric field. We also choose for convenience 
a linear temperature extrapolation for the normal state 
resistivity which vanishes at T = 0, and has a typical 
value p N = 84 (ilJcm at T = 200 K. 

In Figure ^ th e results of the Hartree approximation 
J for the critical fluctuations are compared to the ones ob- 
(43ff^ nec ^ f rom the- Gaussian fluctuation theory. The zero- 
J field critical temperature T c q in the Hartree model was 
considered identical to the mean-field critical tempera- 
ture Tq in the Gaussian approximation, in order to 
have the zero-field transition at the same temperature in 
both theories. This identification causes the mean-field 
transition temperature To in the Hartree model to shift 
upwards with respect to T c o. For our chosen parameters 
this shift was found to be 



T -T c0 = 5.733 K, 



(45) 



To = T c (0) 



2ffT c (0) 



(44) 



In practice, one knows the actual critical temperature 
T c (0) = T c0 measured at very low electrical field, so that 
Eq. I|44|) allows to estimate the bare mean-field charac- 
teristic temperature Tq. Then, having the parameter Tq, 



while taking a cut-off parameter c — 1. The difference 
between the two temperatures depends on choice of the 
cut-off parameter, namely it increases with the c value, 
and becomes divergent for no cut-off (Tq/T c o — ► oo for 
c — > oo). It can be noticed in Figure^that the curves ob- 
tained in the Hartree approximation are less steep than 
those in the Gaussian one, and show a significantly broad- 
ened transition region in the presence of strong applied 
electric fields. In addition, we find that the paracon- 
ductivity in the renormalized model is more sensitive to 
the electric field, showing a more pronounced suppres- 
sion of the fluctuations at high fields in the lower part 
of the transition. However, above the zero-field transi- 
tion temperature, the paraconductivity in the Hartree 
model is always higher than the one in the Gaussian ap- 
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FIG. 1: Resistivity in the Gaussian theory (dotted curves) 
and in the Hartree approximation for the interacting fluctua- 
tions (solid curves), for different values of the applied electric 
field. The following parameters were used: interlayer distance 
s = 1.17 nm; zero-temperature in-plane and out-of-plane co- 
herence lengths, £o = 1-2 nm and £oc = 0.14 nm, respectively; 
Ginzburg-Landau parameter k = 70; zero-field critical tem- 
perature T c o = 92 K. The UV cut-off parameter c = 1 was 
used. The arrow indicates the increasing electric field direc- 
tion. The inset illustrates the critical temperature shift in- 
troduced by the Hartree model if the mean-field transition 
temperature To were kept identical with the one in the Gaus- 



sian approach T 



proximation, due essentially to the critical temperature 
redefinition from Tq to T c q. If one preserved instead the 
same mean-field transition temperature T = T as in 
the Gaussian approximation, one could then visualize the 
critical temperature shift introduced by the Hartree ap- 
proach, as shown in the inset of Fig. ^ Equation 144|) 
would give then T c q — 86.894 K, and the paraconductiv- 
ity would be always lower than the one in the Gaussian 
approximation. 

For illustration, we also give in Figure a comparison 
between the results of our model, applicable for layered 
superconductors in arbitrary electric fields in the absence 
of magnetic field, and the ones of the complementary 
model of Ullah and Dorsey^ which treats in the same 
Hartree approximation the case of an arbitrary magnetic 
field in the linear response (zero electric field) limit. As 
one can easily notice, the same well known "fan shape" 
transition broadening, encountered when a magnetic field 
is applied, can be also predicted for the presence of a 
sufficiently strong electric field. We argue thus that high 
electric fields can be used to suppress order parameter 
fluctuations in HTSC as effectively as a magnetic field. 

This comparison between the fluctuation suppression 
effects of the electric and magnetic field could give also a 
rough estimation on how broad the validity domain of our 
model could be. It is known for instance that the renor- 
malized fluctuation modeU^ for layered superconductors 
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FIG. 2: Comparison between the transition broadening ef- 
fect of the electric (solid lines) and magnetic (dotted lines) 
field, respectively, according to the renormalized Hartree ap- 
proximation. The same parameters as in Fig. 0were used. 
The magnetic field effect was calculated in the linear response 
(zero electric field) limit, according to the UD model. 13 The 
UV cut-off parameter c = 1 in our model corresponds with 
limiting the sum on the Landau levels at the index l/2h in the 
UD model (h — 2iy£qB/&o being the reduced magnetic field). 
The arrow indicates the increasing electric and magnetic field 
direction. 



can successfully fit the resistivity curves of YBCO and 
BSCCO single crystals in magnetic fields up to about 
10 T, in a temperature range that covers approximately 
the two superior thirds of the transition region. This ex- 
tends from a few K in lower fields, up to more than 10 
K in higher fields, measured down from the resistivity 
onset point in zero field^ The lowest third of the transi- 
tion, where, experimentally, the resistivity slope becomes 
steeper, is instead affected by flux pinning effects and 
does not fit into the GL theory. We can therefore as- 
sume, based on the similarities illustrated in Fig. [21 that 
also in high electric fields, the renormalized fluctuation 
model based on the TDGL approach may have its va- 
lidity in a temperature range at least as broad as in the 
case when a magnetic field is applied. Since under high 
electric fields (and consequently, high current densities), 
the pinning of the self-field flux lines is overcome by the 
high Lorentz force, we can expect that the validity of the 
presented model might extend for even lower tempera- 
tures. 

In a few previous papers jSL^S the non-ohmic fluctu- 
ation conductivity in high electric fields was reported 
to be experimentally proven, by confronting the mea- 
sured paraconductivity to the scaling laws predicted by 
Gaussian fluctuation models 4*^ However, the broaden- 
ing of the temperature dependence of the resistive tran- 
sition with respect to the increasing electric field, and 
the breaking of the mean-field (Gaussian) theory in the 
immediate vicinity of T c , signalized by Ref. l39i indicate 
that a renormalized (non-Gaussian) fluctuation model, 
as the one presented in this paper, might be more ap- 
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propriate. From the experimental viewpoint, applying 
electric fields of a few hundreds V/cm on cuprates may 
however be not an easy task, since the dissipated power 
density would attain levels of the order of GWcm~ 3 . On 
the one hand, high electric fields are necessary in order to 
put into evidence the non-ohmic fluctuation conductivity, 
while, on the other hand, they produce high dissipation 
and can increase the sample temperature at values where 
the nonlinearity is no longer discernable. In this connec- 
tion, using short current pulses at high current densities 
(a few MAcm~ 2 ), seems to be a better alternative to the 
dc and ac measurements. 



more pronounced suppression of the fluctuations at high 
fields in the lower part of the transition, and a higher 
paraconductivity above the zero-field transition temper- 
ature than the Gaussian fluctuation model. All these 
features are quantitatively important for commonly used 
HTSC, so that experimental investigations could be able 
to discern easily between the applicability of this model 
in competition with the Gaussian fluctuation approxima- 
tion. 
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VI. CONCLUSIONS 

In summary, we have treated in this paper the critical 
fluctuation conductivity for a layered superconductor in 
zero magnetic field, in the frame of the self-consistent 
Hartree approximation, for an arbitrary electric field 
magnitude. The main results of our work are the formu- 
lae ()24[) for the fluctuation conductivity, and (|29|) for the 
renormalized reduced-temperature parameter. In both 
equations the UV cut-off of the Ginzburg-Landau model 
was explicitly considered. In the linear-response limit 
(E — ► 0), the corresponding expressions Eqs. 131|) and 
(|S2^1 reduce to the previous results of existing theories. 
Qualitatively, the temperature characteristics at different 
electrical fields in the Hartree approximation turn out to 
be less steep than those in the Gaussian one, they show a 
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APPENDIX: GREEN FUNCTION FOR THE 
TDGL EQUATION 

Equation (|12f) can be solved easier for the Fourier 
transform of the Green function with respect to time, 



R q {k,uj;k' x ,t') = JdtR q {k,t; k' x ,t') 
which satisfies the equation 



(A.l) 



«i 



h dk x 2m 
One finds for the differential Eq. (|A.2|I the solution 

R q (k,Lo;k' x ,t') = A q (k,uj;k x ,t')-exp 



R q (k,cj;k' x ,t') = 6(k x -k' x ) 



\e E 



— h (ai - iujT ) k x 

6m 



where the derivative of the coefficient A q (k, u; k' x , t') must satisfy 



dk x 





exp< 


r hr 


\h 2 kl 


e E 


{ e E 


6m 


nr 


exp< 


r ^r 


'h 2 k' x 3 


e E 


I e E 


6m 



— + (ai — iujT 1 ) k a 



+ (ai — iuiT 1 ) k' x 



5{k x - k' x ) 
S(k x - k' x ). 



The solution for the coefficient A q (k,Lu; k' x ,t'), which remains nondivergent when k' x — > — oo, is 



A 9 (k, W ; k' x ) = ^-9{k' x - k x ) exp I 



hk' x 3 
6m 



fax — iujY 



J n x 



e E y x r I e E 

where 6(k' x — k x ) is the Heavyside step function, so that the Fourier transform of the Green function will be 

^ 2 ~ k'x) 



R q ( k ,^k' x ,t') = ^e(k' x -k x )-e X pl^ 



6m 



+ (a x - iuT Q : ) (k x - k' x ) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 
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Now we can apply the inverse Fourier transform to regain the Green function depending on time i, and obtain 



R q (k,t; k x ,t) 



2vr 



e E 
2^ 



v (ki-*g) 



6m 



ai (k x - k' x ) 



exp < — no 



t — t-\ — {k x — k ' ) 
e E 



= -^-^ 6(k' x - k x ) exp \ — ^ 



e E 



e E 



h 2 (k^. fc^ 3 ) 



6 m 



+ ai (k x - k' x ) 
I 



(A.7) 



S[t-t' + — [k x -k' x ] 



The form (|A.5|) for the coefficient A q (k, to; k x ), and (|A.6|1 doesn't diverge, and provides also the retarded 
namely the presence of the Heavyside function 6(k' x — k x ) character in Eq. l|A.7jl . i.e. R q (k,t; k' x , t') = for t < t'. 
assures that the Green function i? g (k, lo; k' x , t 1 ) in Eq. 
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